Self-accelerating Massive Gravity: 
Time for Field Fluctuations 
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The ghost-free theory of massive gravity has exact solutions where the effective stress energy gen- 
erated by the graviton mass term is a cosmological constant for any isotropic metric. Since they are 
exact, these solutions mimic a cosmological constant in the presence of any matter-induced isotropic 
metric perturbation. In the Stiickelberg formulation, this stress energy is carried entirely by the 
spatial Stiickelberg field. We show that any stress energy carried by fluctuations in the spatial field 
away from the exact solution always decays away in an expanding universe. However, the dynamics 
of the spatial Stiickelberg field perturbation depend on the background temporal Stiickelberg field, 
which is equivalent to the unitary gauge time coordinate. This dependence resolves an apparent 
conflict in the existing literature by showing that there is a special unitary time choice for which the 
field dynamics and energy density perturbations vanish identically. In general, the isotropic system 
has a single dynamical degree of freedom requiring two sets of initial data; however, only one of these 
initial data choices will affect the observable metric. Finally, we construct cosmological solutions 
with a well-defined perturbative initial value formulation and comment on alternate solutions that 
evolve to singularities. 
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I. INTRODUCTION 

Constructed to remove the Boulware-Deser ghost [l|, 
the theory of massive gravity also possesses solu- 

tions that accelerate the cosmological expansion in the 
absence of a true cosmological constant [a-flg . In a pre- 
vious paper (l3j ]. we demonstrated that all of these ap- 
parently distinct solutions are part of a single class of 
isotropic solutions in the Stiickelberg formulation of the 
theory. This class possesses several key features. The 
spatial Stiickelberg fields follow the radial coordinate of 
the metric itself and produce the same effective cosmolog- 
ical constant in the presences of any isotropic metric. The 
temporal Stiickelberg field is inhomogeneous in isotropic 
coordinates but does not play a direct role in establish- 
ing the effective cosmological constant. Solutions for this 
field are not uniquely specified, leading to seemingly dis- 
tinct versions of self-acceleration. 

Once the Stiickelberg fields arc set on a self- 
accelerating solution, they remain on it for any evolu- 
tion of the matter fields that remains isotropic, includ- 
ing cosmological expansion and radial collapse of mat- 
ter perturbations. Given the inherent interest of self- 
accelerating cosmological solutions, it is important to un- 
derstand whether Stiickelberg field perturbations around 
them are themselves stable and healthy. This question 
has now been widely studied, with different authors us- 
ing different approaches. Confusingly, these different ap- 
proaches seem to give different conclusions. Study of 
perturbations in a decoupling limit showed a potentially 
healthy scalar degree of freedom but potentially prob- 
lematic (strongly coupled or ghost-like) vector degrees of 
freedom 0, [lji [l8| ■ A local patch expansion approach, 
which one might expect to be qualitatively similar to a 
decoupling limit, found no propagating Stiickelberg de- 
grees of freedom at all [l9[. Finally, study of perturba- 



tions in the full theory around a particular solution for 
the temporal Stiickelberg field similarly found no propa- 
gating degrees of freedom (2(| ■ 

To resolve this issue, we study spherically symmetric 
Stiickelberg perturbations around the full theory for the 
whole class of solutions. We find that in general the 
spatial Stiickelberg field fluctuation does possess dynam- 
ics and carries stress-energy fluctuations in addition to 
the background constant. There are special choices of 
the temporal Stiickelberg field background that elimi- 
nate the dynamics but also eliminate the possibility of 
energy density fluctuations. In all cases, isotropic energy 
density fluctuations are stable and from general initial 
conditions decay back to the cosmological constant in an 
expanding universe. Our approach allows us to show that 
the decoupling limit findings [(| [ItJ and specific temporal 
solutions in the full theory [20j | are not, after all, in con- 
flict. We have not been able to harmonize our findings 
with Ref. fl9j , who finds no dynamics for either temporal 
Stiickelberg field choice. 

The structure of this paper is as follows. In Jl] we 
briefly review the theory of massive gravity [5|| a nd the 
class of isotropic self-accelerating solutions [13} • We 
then consider the equations of motion and action gov- 
erning spherically symmetric field perturbations around 
this class of solutions in mill Finally in SjTV] we show 
how the choice of the temporal Stiickelberg field in the 
background affects the dynamics of the perturbations and 
their contribution to stress energy. We close with a dis- 
cussion and summary of these results in ifVl 



II. SELF-ACCELERATION 

In the section, we review the action, equations of mo- 
tion, and stress energy for the massive gravity model 
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[H HI • We specialize these quantities for the class of exact 
isotropic self-accelerating solutions [l3[ which will form 
the basis of the perturbation studies that follow. 

A. Massive Stress Energy 



The Lagrangian density [H, H| 



Ml 



(i) 



represents a covariant theory of massive gravity con- 
structed so as to eliminate the Boulware-Descr ghost 
P, Hi], Hit ■ Here m is the graviton mass and M v \ is the 
reduced Planck mass. The potential IA is constructed in a 
covariant manner by making it a function of the so-called 
fiducial metric E, 



S, 



(2) 



where 4> a are the 4 Stiickelbcrg fields. These fields trans- 
form as spacetime scalars and hence maintain general 
covariance. 

The potential hi can now be written in terms of the 
matrix 



^ V — 9 Z-'OLV, 



(3) 



as 



j = -12 + 6[\/S] + [S]-[\/S] 2 
+ a 3 (-24+18[VS] -6[Vi:}'- 
-3[S]([v / S]-2) + 2[S 3 / 2 ] 

+ a 4 (-24 + 24[v / £] - 12[%/S] 2 - 12[VS][X!] 
+ 6[%/S] 2 [S] + 4[\/S] 3 + 12[S] - 3[£] 2 
- 8[S 3/2 ]([7S] - 1) + 6[S 2 ] - [x/S] 4 ) , (4) 



where brackets denote traces, [A] = A* 1 ^, and 03, a 4 are 
free parameters. 

Variation of the action with respect to the metric yields 
the modified Einstein equations 



m T„ 



_}_rp{m) 

Af 2 ^ ' 



(5) 



where G^ v is the usual Einstein tensor and is the 

matter stress energy tensor. Here 



1 



.U 



-9 Sg» v 



(6) 



is the dimcnsionless effective stress energy tensor pro- 
vided by the mass term. An explicit expression in terms 
of X is given in Ref . [l3[ . 



B. Effective Cosmological Constant 

A constant stress energy is an exact solution of massive 
gravity for any spatially isotropic metric fl3j , 

ds 2 = -b 2 (r, t)dt 2 + a 2 {r, t)(dr 2 + rW). (7) 

To see this fact, consider that an isotropic parameteriza- 
tion of the Stuckelberg fields 



0° = /(t,r), 
! =g(t,r) — , 



(8) 



enables the potential to be written more compactly as 



\ar 



WB: 



(■*■), (9) 
\arJ 



where 



MM) 2 -(f 



A 4 — , 
a 



(10) 



and /x = sgn(/</ — 5/'). Here and below overdots denote 
derivatives with respect to t and primes denote deriva- 
tives with respect to r when acting on Stuckelberg or 
metric fields. Note that W is related to the determinant 
of S 1 / 2 which we assume is never zero. Correspondingly 
fi must be either +1 or —1 throughout the spacetime. 
This will be an important consideration in QYV\ 
The P n polynomials are 

P (x) = -12 - 2x(x - 6) - 12(x - l)(a; - 2)a 3 

- 24(x- l) 2 a 4 , 

Pi (a;) = 2(3 - 2x) + 6(x ~ l)(x - 3)a 3 + 24(a; - l) 2 a 4 , 

P 2 (x) = -2 + 12(x - l)a 3 - 24(x - l) 2 a 4 , (11) 

and satisfy the recursion 

P' n =2P n+l -xP' n+1 , (12) 

where here and throughout P^(x) = dP n /dx and should 
not be confused with radial derivatives. 

Varying the action with respect to / and g yields the 
Stuckelberg field equations 



0, 



X ^ a ■ 



and 



— d r 



-di 



a 2 br 2 
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X \b 
a 2 br 2 



L 

a 



\ia 2 r 2 g'P 2 



Pi + fia 2 r 2 gP 2 



(13) 



= 0, 



fi J -)P 1 + f xa 2 r 2 f'P 2 
a ' 



Mt + - Pi 

b a . 



Ha 2 r 2 fP 2 



= a 2 br \P^ + VXP{ + WP2 



(14) 
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Eq. (fT3|) has an exact solution given by Pi(xq) = 0. Here 



1 + 6c*3 + 12a 4 ±Jl + 3a 3 + 9a\ - 12a 4 . , 
*° = 3 (a3 +4a 4 ) (15) 

unless a3 = a4 = in which case xq = 3/2. Given 
g = xgar, Eq. (|14p reduces to 



HP W 
V A = h Xo, 

a; 



(16) 



and can be used to define solutions for / (see glVJ) . 

For any such solution, the stress energy tensor takes 
the cosmological constant form 



where 



-P = 2 P o(^o)- 



(17) 



(18) 



To restore physical units to the dimensionlcss stress ten- 
sor multiply by m 2 M^. The Einstein equations (|5|) then 
determine the metric functions a and b jointly in the pres- 
ence of matter in a manner independent of the solution 
for /. The self-accelerating stress energy depends only 
on the universal form for the spatial Stiickelberg field, 
g = XQar. common to the whole class of isotropic solu- 
tions. 



III. FIELD FLUCTUATIONS 

In this section we derive a complete set of equations of 
motion for the observable properties of Stiickelberg field 
fluctuations or equivalently the field equation for spatial 
Stiickelberg fluctuations. We then relate these results 
to the second order action and its dependence on the 
temporal Stiickelberg background solution and massive 
gravity parameters. 



A. Equations of Motion 

Now let us consider spherically symmetric Stiickelberg 
perturbations. Hereafter / and g represent the back- 
ground solution for the Stiickelberg fields, a and b the 
background metric, while Sf and 



5T = 5g — xorSa, 



(19) 



5a and 5b quantify the perturbations. Because our so- 
lutions are exact for any a(r, t), the metric can have 
perturbations away from a given background, e.g. the 
Friedmann-Robertson- Walker background, due to per- 
turbations in the matter sector and still be on the ex- 
act self-accelerating solution in the Stiickelberg sector. 
5T therefore is the fluctuation in g away from the self- 
accelerating solution. 



Linearizing Eq. (|13|) in the fluctuations, we obtain a 
closed-form equation for ST 




— d r 



abr 

7x 



t 

a 



ST 



(20) 



This equation for ST contains no reference to matter 
or metric perturbations. It is also decoupled from the 
Stiickelberg fluctuation Sf. The coefficients in this equa- 
tion do however depend explicitly on the background so- 
lution for the metric and both Stiickelberg fields. 

This equation is first order in time and requires one set 
of initial data ST(r, t = 0) to solve. However, note that if 
we choose initial values of ST(r, 0) = 0, ST will stay zero 
regardless of matter-driven perturbations in the metric. 
In this sense, the Stiickelberg fluctuations are decoupled 
from the matter. 

As was the case for the background, the perturbed 
Eq. (fT4"| defines the evolution of Sf and depends explic- 
itly on ST as well as on the metric fluctuations (see Ap- 
pcndix[3]) . Thus to characterize the complete Stiickelberg 
sector requires a second set of initial data 6f(r,0). Con- 
sidered together, the amount of initial data and the two 
coupled first order systems are indicative of a single prop- 
agating degree of freedom. However, these fields are not 
interdependent in the usual way of a field and its conju- 
gate momentum and in particular do not combine to form 
a single wave equation: ST does not depend on Sf. We 
will return to the interpretation of these facts in 3111 Bl 

Even though the Stiickelberg fluctuations are decou- 
pled from the matter, a finite ST generates a metric per- 
turbation, and hence a matter perturbation, due to the 
effective stress energy it carries. By explicitly evaluating 
Eq. ([6]), we can write down the energy density perturba- 
tion 

~' 2 - f' 2 + a 2 W _ \ P[5T 



Sp = -ST° = 



Xq 



IX 'J 2ar 

and the radial momentum density perturbation 



Sq = T' 



ff - 99' PjST 



■, (21) 



(22) 



Here and below the polynomials Pi = P,(xo) are assumed 
to be evaluated on the background. Note that both quan- 
tities are directly proportional to P[. Furthermore Sp and 
Sq depend only on ST and not its derivatives, in contrast 
to the energy and momentum of usual fields, which typi- 
cally also get contributions from the their kinetic energy 
The energy and momentum density suffice to define the 
impact of the Stiickelberg fields on matter through the 
two associated Einstein equations. To define the whole 
stress tensor or use all of the Einstein equations, we can 
also derive the radial pressure 



ST r 



r- 



b 2 W 



b 2 Vx 



Xq 



P[ST 
2ar ' 



(23) 
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The other components of the stress energy tensor can 
then be found through conservation of energy 



(6a 



3V 



6 P = -3-(6 P + 8p) + ^ 



Sq- 



(r 2 5q)' 



0. (24) 



Given Sp and Sq, this equation defines Sp. Since by defi- 
nition 



5p=\{5T\ + 5T% + 5T\) 



(25) 



and ST 9 9 = ST^ , the energy equation defines the remain- 
ing components. We have verified that direct evaluation 
of the stress components in Eq. (J6j) gives the same result 
once combined with the equations of motion. With the 
stress tensor fully defined, any two of the Einstein equa- 
tions completes the dynamics as usual and define 5a and 
8b jointly with the matter stress energy. 

Momentum conservation provides another useful aux- 
iliary equation 



Sq 



/a 6\ 6 2 
5---Ug+ — 
V a 6/ a 1 



h -[t 



(26) 



-{ST r r y 



(ar)' 



28T\ - 5T e n - STt 



With an anisotropic stress tensor, balancing radial gradi- 
ents with the different components can result in static 
conditions if 



(8T r r 



(ar)' 



ar 



25T r l - ST e a - ST 4 



(27) 



Whereas for an isotropic stress tensor, the right hand 
side vanishes and (5T r r )' = Sp' , whose finite value would 
generate a momentum density unless balanced by the 6'/6 
gravitational term. The impact of anisotropic equations 
of state on hydrostatic equilibrium and stability has been 
studied in a related neutron star context (e.g. [23l [24l|). 

In summary, similar to what we found for the back- 
ground solution, the stress energy tensor for the fluc- 
tuations depends on the perturbation to the spatial 
Stiickclbcrg field (<5r) and not the perturbation to the 
temporal Stiickclbcrg field (5f). Hence even though the 
Stiickclbcrg fields require two sets of initial data obeying 
coupled first order equations, reminiscent of a field and 
its conjugate momentum, only one of these has any ob- 
servable impact on the matter. Moreover, the coupling is 
unidirectional, as ST forms its own autonomous first or- 
der system. It can be consistently set to zero by a choice 
of initial condition. 



B. Action 

In order to make contact with the literature on 
Stiickclbcrg field dynamics, it is useful to see how these 
results arise from the second order action. In particular, 
apparent differences concerning the sign and value of the 



time derivative terms in this action have led to seemingly 
conflicting results about the presence of strongly coupled 
or ghost- like modes @, [l?], Ho] ■ 

We expand the action of Eq. © to second order in 
both Stiickelberg and metric perturbations around our 
background solution in the Appendix. It is sufficient here 
to examine the Stiickelberg sector: 



Sss 



drd£L 



SfST eom - -ST 2 



2 br8T(8T'd g > +STdg)VX 



(28) 



where oT oom was defined in Eq. ([2H| and note that this 
is the only term that depends on 5f in the whole sec- 
ond order action (see Eq. IA5|) . Furthermore since <5r eom 
depends only on <5r, there are no terms quadratic in Sf 
fluctuations. Varying the action with respect to Sf re- 
produces the perturbed equation of motion $r oom = 
with no dependence on other perturbations. 

With the help of Eq. (fTB")) , we can see that the coeffi- 
cient of the SfST term in the action is proportional to 



bW 



(ar)> 



(29) 



and determines the dynamics of ST through its equation 
of motion. The quantity in brackets is determined by 
the background solution for / and there is a special case 
where it vanishes identically. Note that if we assume 
that the time dependence of / cx a p , then / oc fa/ a and 
bW cx fa/ a so this condition becomes independent of 
time. We shall see that this very special f-g symmetry 
is exploited in open universe solutions (sec ^IVC|) and 
explains their lack of Stiickelberg dynamics |20j . 

While P[ in Eq. (f2T)|) drops out of the equation of mo- 
tion, it does determine the sign of all Stiickelberg coeffi- 
cients in the action. For the self-accelerating solutions 



P{(x )=±4jl + 3a 3 



9o§ 



12a 4 , 



(30) 



for the two Xq solutions. If a 3 = = there is only one 
solution, P[ = —4 and the sign of the coefficient is fixed 
entirely by the background solution. There is a special 
choice 



04 



1 

12 



(l + 3a 3 + 9a§), 



(31) 



where P[ = 0. In this case the whole Stiickclbcrg 
quadratic action vanishes as do all of the stress-energy 
components. This is the same special choice that was 
made in Refs. [1, [l(|, who similarly found that the 
quadratic action for fluctuations around solutions with 
this parameter choice vanished. 

For general a 3 and 0:4, P[ and hence the action may 
have either sign. This is consistent with the decoupling 
limit analysis of the kinetic term where a 3 7^ or a± 7^ 
was necessary to ensure that the helicity scalar fluctua- 
tion is not a ghost [H,[l7]]- In this limit, the scalar, Sir, de- 
termines the two Stiickelberg field fluctuations ST — » Sir' , 
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Sf —> —Sir and hence the SfST term acts as a kinetic term 
for Sir after integration by parts (see also Eq. l4"Tj) . 

The decoupling limit form suggests a partial way to 
interpret the first order nature- of the coupled Sf-ST sys- 
tem. There Sf appear in the action in a way similar to 
a canonical momentum for the field Sir, which is itself 
related to ST [25| . This viewpoint suggests an interpre- 
tation of Sf and ST as a pair of "half degrees of freedom. 
This counting is related to the absence of the Boulwarc- 
Deser ghost: of our original 4 Stiickelberg fields, there 
is one constraint by construction and there are two an- 
gular modes that are absent in our analysis because of 
spherical symmetry. 

However, interpreting Sf and ST as a simple single 
propagating scalar does not carry over to the full the- 
ory, where the fields are not directly related. This is of 
course not in itself a contradiction, as the structure of 
this theory makes a global helicity decomposition gener- 
ically impossible [2tSj . In particular Sf has no necessary 
relationship to ST as any choice that solves its equation 
of motion yields the same solution for ST. Furthermore, 
the dynamics of ST depends on the choice of background 
solution for /, as we shall see in the next section. 



IV. TIME DYNAMICS 

The detailed dynamics of Stiickelberg field fluctuations 
depend on the choice of solution for /, the Stiickelberg 
field related to the time coordinate. There are even spe- 
cial cases where the fluctuations have no dynamics. All 
such choices have the same background metric and self- 
acceleration. In this section we systematically construct 
explicit solutions for / and study their implications for 
the dynamics of field fluctuations and the associated ini- 
tial value problem. We also generalize many solutions 
in the literature HH, [lTj as well as introduce new 
classes. 



A. Unitary Gauges 

Since the Stiickelberg fields are spacetime scalars, their 
values, once mapped to the same spacetime point, are 
coordinate independent. Unitary gauges arc defined so 
their time t u and radial r u coordinates are 



tu = f{r,t), 



g(r, t). 



(32) 



The Stiickelberg fields are therefore just the unitary 
gauge coordinates and their functional form in isotropic 
coordinates (r, t) is simply the map itself. In unitary 
gauges, the fiducial metric S M „ = -q^ and this simplifi- 
cation has been crucial for finding most solutions existing 
in the literature. While we do not utilize a unitary gauge 
construction in our solution, we will categorize solutions 
by their unitary gauge correspondence. 

All self-accelerating solutions have the common prop- 
erty that g = XQar, implying that the radial coordi- 



nate of any unitary gauge is simply the same conformally 
rescaled version of the isotropic radial coordinate, similar 
to the distinction between physical and comoving coor- 
dinates in cosmology. 

The different solutions are thus characterized by the 
choice of time coordinate for the unitary gauge. Each 
of these choices solves Eq. (fl6|) for the same g; written 
explicitly, 



b 2 f' 2 + 2ar(a'f 2 aff) + r 2 (a'f - aff = 
x 2 (a' 2 b 2 r 2 + 2a'ab 2 r-d 2 a 2 r 2 ). 



(33) 



Recall that the metric functions a and b do not depend 
on the solution for / and so may be considered external 
functions in solving this equation. 

There are several properties of this equations that are 
useful to note. It is a nonlinear partial differential equa- 
tion in space and time. Its solutions can be characterized 
by the boundary conditions /(0, t) and f(r, 0). The for- 
mer is equivalent to the local r = relationship between 
t and unitary time / and is our primary classification 
criterion. The latter allows a family of such solutions as- 
sociated with a time integration constant. The equation 
does not guarantee that the map in Eq. (|32|) is free from 
singularities everywhere in the spacetime. A singularity 
corresponds to a change in the sign of the Jacobian de- 
terminant fg' — gf and hence a change in the sign of 
/i (see Eq. [T0|) . Solutions should be checked against sin- 
gularities developed in the course of radial or temporal 
integration as they signal a breakdown in our treatment 
of the square roots in the action. 

Finally given a solution for the = a.4 = case where 
x = 3/2, 



f(r,t;ot3,an) 



2x 



/M;0,0) 



(34) 



is the solution for the full parameter space. One may 
similarly scale any other specific 0:3 and on case to the 
general case by use of the ratios of the respective values 
of xq. This relationship makes it trivial to generalize 
special solutions in the literature (cf. @). 

In our examples, we seek cosmological solutions where 



b(r,t) = l, 



a(r,t) 



a F (t) 



1 + Kr 2 /4' 



(35) 



where K is the spatial curvature and ap(t) is the scale 
factor of a Friedmann-Robertson- Walker metric. Note 
that in a flat K = geometry a(r, t) = a(t) = a^(i). 



B. Unitary Time t 

Perhaps the most natural choice of solution is to take 
unitary time / to be linearly related to t 



f(0,t) 



£0 
C 



t. 



(36) 
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where C is a constant. For a flat universe K = and we 
can solve Eq. ([55)) with this ansatz to 0{r A ) 



f(r,t) = ^lm + 1 -^^l (aHrr 



(37) 



where H = a /a. Note that on cosmological scales where 
aHr > 1, this temporal field becomes spatially inhomo- 
geneous. 

1. de Sitter 

The special case that H = const, is of course the late 
time limit of the self-accelerating solution where the mat- 
ter becomes subdominant to the effective cosmological 
constant 



W-Pom, 



(38) 



and we can normalize a = e Ht . Solving Eq. (|33j) for this 
case with the ansatz 



f(r,t) = -^[Ht + F(aHr)\ 



(39) 



gives 



Ht-y + In 



1 + y 



1 - {aHr) 2 



y = ^Jl + C 2 {a 2 H 2 r 2 - 1). 



(40) 



These relations are equivalent to the de Sitter solutions 
of Refs. @, HI up to an overall constant and generalizes 
them to arbitrary a^, a^. 

Note that if C = 1, / oc t and g oc ar = r p near r = 
and can be thought of deriving from a Lorentz scalar in 
physical coordinates 



1 



(41) 



with / w — 7r and g sa <9r p 7!"- Hence C ^ 1 may be said 
to have a "vector" in the background || even though the 
Stiickelberg fields are spacetime scalars. 

Given the solution in Eq. (|40p . the evolution equation 
for <5r becomes 



(y 2 + y) 



ST 
~H 



y 



{aHr) 2 

2y 2 



rST' 



[aHrf 



iy l 



l 



sr. 



(42) 



and the energy density, momentum density and pressure 
are given by 



6p = - 



(aHr) 2 C 3 x P[5T 
(1 + C)(1 -y) 2 2ar ' 



Sq 



y + c 2 

C 2 a 2 Hr 



1 



-Sp, Sp = —Sp . 



(43) 



The pressure follows a radiation-like equation of state but 
itself is composed of anisotropic contributions 



y + c 2 



1 



C 2 aHr 



Sp 7^ Sp. 



(44) 



It is instructive to consider the simplest C = 1 case 
further. Here y = aHr, and given an initial value for 
ST(r, 0), the general solution is 



c , . If a-l\ 3 / a- 1 

which then defines the stress components 
{aHr) 2 x P{ST 



(45) 



2(1 - aHr) 2 2ar 



1 



<5g = —5/0, (5T r r = 'iSp = Sp. 



(46) 



The angular stresses vanish and the total pressure is 
given by the radial component. These components have 
poles in their expressions at the horizon r = 1/ aH corre- 
sponding to the coordinate singularity in Eq. (|40|) for 
unitary time. However, since the solution for ST al- 
ways maps the initial horizon onto the horizon at a later 
epoch, if the energy density is finite in the initial condi- 
tions there, it remains so. Furthermore for large scales 
r 1/H, which were outside the horizon at the initial 
epoch, Sp oc a -4 ; this is consistent with expansion effects 
and the radiation-like equation of state. 



2. Matter and Radiation Domination 

While we have solved the initial value problem for fluc- 
tuations in the de Sitter limit for f(0,t) oc t, in a cos- 
mological solution these would themselves originate from 
the dynamics in the preceding radiation and matter dom- 
inated epochs. 

To construct solutions where /(0, t) oc t during the 
matter and radiation epochs, we first try 



f(r,t) = ^F(ar/t), 



(47) 



with the assumption that H oc a~ 3 ( 1+w ^ 2 for a constant 
w and C. For any such choice, we can solve differential 
equation for F implied by Eq. ([55| . While this construc- 
tion would seem to admit many solutions, for finite C, p, 
changes sign between +1 for r — > and —1 for r — > oo 
indicating a singularity in the map between isotropic and 
unitary coordinates. 

On the other hand, the C — > oo limit does not suffer 
this problem and yields 



f{r,t)=x ^/r 2 {t)+r 2 a{ty 



where 



/ 9(l + w) 2 
T ~ V 4(2 + 3w) 



(48) 



(49) 
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for w > —2/3. Eq. (|49[) is now easy to generalize to an 
arbitrary expansion history by matching 



(50) 



for d > 0. This solution includes a universe that evolves 
from radiation domination through matter domination 
to self-acceleration. Note that since w = — 1 during self- 
acceleration, t is no longer directly related to t and so 
this solution is distinct from the solutions of the previous 
section. 

The fluctuation ST evolves under 



2rVd<5r + [dV 



, 6V > = -2pr, (5i) 



where 



A = d 4 r 4 + 2ad J < + rVdV - 2a J dT r (l - r 2 d 2 ) 

+2r 2 a 5 d - a 4 (l - dV + 2d 4 r 4 - 2dd> 2 7>), 
B = (a + dr r ) 2 - rVd 2 , (52) 

with t 2 == t 2 -(- r 2 a 2 and the stress energy components 
are 



a 2 dr xoP^r 



Sq 



2ra 3 d 



-<5p, 



£p = ^ r 



2 ■ 2 2 \ 2 

a — arr 



2ra 2 c 



(53) 



Note that for a constant background equation of state 
w, 5p/Sp = —(1 + 3w)/6. For w > — 1/9 expansion thus 
makes the energy density redshift more slowly than the 
dominant matter but always faster than the constant 
background Stiickclbcrg contributions. Thus in linear 
theory starting from some arbitrary field configuration, 
we expect that Sp/p will be driven rapidly to zero for 
any matter content leaving just the cosmological constant 
background. 



1. Open Solution 

If K < then Eq. (|55[) has a simple solution for C 2 
-K. Taking this ansatz, the full solution is [l2l. [20j 



f(r,t) = x a F (t) 




(1 + J0 2 /4) 2 



(56) 



Note that the open solution for /, like other solutions, is 
inhomogcncous in isotropic coordinates. There is nothing 
special about the open solution with regards to homoge- 
neous and isotropic Stiickelberg fields (cf. fl2l]). 

On the other hand, the common separable ap{t) factor 
in /, g and a allows this solution to satisfy the static 
condition for the <5r field, Eq. (|29p . Correspondingly, 
Eq. ([20]) becomes 



(4 + Kr 2 )r5Y' + 2(4 - Kr 2 )5Y = 0. 



(57) 



Note that p = sgn(dF) and so is +1 for an expanding 
universe. Interestingly the determinant goes to zero if 
the expansion turns around, e.g. because of a negative 
true cosmological constant, signaling a breakdown of the 
solution. We therefore consider only expanding solutions 
here. 

The general solution to this equation is 



ST(r) ex 



(58) 



This static Stiickelberg field produces no energy density, 
momentum density or pressure; i.e., 



Sp = Sq = Sp = 0. 



(59) 



These static conditions are maintained by a delicate bal- 
ance of the radial and anisotropic stress gradients 



ST r 



1 



hp \ x$P[5T 



2ar 



(60) 



C. Unitary Time a 

Unitary time can also be made proportional to 
a(Q,t) = a F {t) 



x 



f{0,t) = -ga F (t) 



(54) 



where, again, C is a constant. In this case both / and 
g share the same temporal dependence at the origin and 
hence provide the ingredients necessary for static solu- 
tions. To order 0(r 4 ) 



(55) 

Here we have kept the possibility that K ^ since it 
allows a special class of solutions. 



which satisfy the static condition for Eq. ([27|) . 

Even though Sp = Sq = 0, anisotropic stress still has 
an impact on the metric through the Einstein equations. 
Moreover, these types of solutions are potentially unsta- 
ble to anisotropic perturbations [13, HH . 



2. Flat Solution 



UK = 0, then there is still one simple solution to 
Eq. ((55)) where C 2 <C a 2 F . In that case 



/(r,t)«^a(t). 



(61) 



To promote this to an exact solution we follow Ref. [ll| 
and generalize their approach to arbitrary 03 , 0:4 



1 



C 2 2 C 2 T 2 {t) 



4 a 2 (t) 



(62) 
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which recovers the approximate scaling for r< 1/C and 
t < a/C. Eq. (HOI) becomes 



2C 2 r 2 aaST + [C 2 + (C*V - A)a 2 ]r8T' = 2— ST (63) 



B 



with 



A = pC 4 + 4C 3 d - I6Cd 3 - p{16 - C 4 r 4 )d 4 



-2CVa(C/i + 2d)a, 

-/xC* 2 - 4C*d - /x(4 - C 2 r 2 )d 2 



(64) 



Despite solving Eq. (f3"3")l . there is a problem with this 
form as a global solution since 



/' 



1 aHr < Jl + AH 2 /C 2 



-1 aHr > Jl + AH 2 /C 2 



(65) 



The sign change in p indicates the mapping to unitary 
gauge is singular unless C — > 0. 

In this limit we are again left with a static equation 



rSV = 28T, 



(66) 



whose solution is <5r oc 1/r 2 . Thus the K — > flat limit 
of the open solution is the same as the C ~ > solution 
of the flat solution. More generally, the stress-energy 
components become 



CVa x P[ST 

°P = n «77 J d ° 



C 2 -(A-C 2 r 2 )a 



D 



2ar 



2C 2 ? 



■Sp, 



C 2 - (4-C 2 r 2 )d 2 



2C 2 rd 



8p. (67) 



Again for the C — > case, <5p, i5g, <5p are all suppressed 
whereas ST r r , ST e e , ST^ are unsuppressed and contain 
radial gradients that delicately balance each other. 



D. Unitary Time y/a 

While /(0, t) ex t or a are obvious choices for solutions, 
they do not exhaust the possibilities. We have seen that 
the former generally require dynamical Stiickelberg fluc- 
tuations whereas the latter admit static ones. However 
those static cases are only formally well-defined as global 
solutions for K < 0. 

To see whether there exists a well-defined static solu- 
tion for K = we construct solutions where /(0, t) oc y/a 
is the limiting form. Similarly to the flat ex a case, we 
find an exact flat solution 



f(t, r) = x 



a(t) 
C 2 



l {t)+a 2 (t)r 2 



(68) 



for any constant C. Note that if C 2 3> r 2 /a, this solution 
recovers the general /(0, t) oc t solution of i jlVB 21 Since 
we have already considered the general case of this limit, 
it is interesting to consider the opposite one. Note that 
in a cosmological solution, once self acceleration sets in 



t 2 I a — > const, of order H~ 2 . Thus taking C <C m ensures 
that this limit is satisfied for all time. 

In this Ccm and a = e Ht case, to leading order 5T 
obeys 



2ar 2 C 2 



ST 



rST' 



-25T, 



(69) 



which is a stiff equation as r — > with an equilibrium 
static solution of ST oc 1/r 2 . The stress components 



5p = - 



2 C 3 x P{5T 
~H Jar 



S P =-5p, 5T\ 



Sq 

H 2 
Ar 2 C 4 



H 



2arC 2 



5p, 



Sp. 



(70) 



Since C <C H, ST r r ^> 6p and again involves a delicate 
balance of anisotropic stresses. This case is very similar 
to the C — s- case of the previous section but has the 
benefit that p = 1 everywhere and so C can be set to a 
finite number. 



V. DISCUSSION 

We have presented a general analysis of isotropic 
Stiickelberg field perturbations around the full class of 
self-accelerating solutions [l3| of the massive gravity the- 
ory. These background solutions arc defined by two 
fields, one spatial g(r,t) and one temporal f(r, i), where 
only the spatial one is responsible for the stress-energy 
of the effective cosmological constant. 

Likewise Stiickelberg field perturbations come in two 
classes, spatial (<5r) and temporal (5f). Spatial pertur- 
bations can be consistently set to zero by a choice of ini- 
tial conditions. They are not generated by any matter- 
induced metric perturbations as was already apparent 
in that the self-accelerating solution is exact and non- 
perturbative for isotropic metrics [l3j . With an arbitrary 
choice of initial conditions, spatial Stiickelberg fluctua- 
tions generically possess stress energy. This stress en- 
ergy produces metric fluctuations to which the matter 
responds. Temporal Stiickelberg fluctuations carry no 
stress energy and have no effect on metric fluctuations at 
this order in perturbation theory. 

Importantly, the dynamics of spatial Stiickelberg fluc- 
tuations make energy density deviations from the con- 
stant background always decay with the expansion in 
cosmological solutions, implying that the background so- 
lution is stable. There are special choices where the 
Stiickelberg fields perturbations arc static, but in those 
cases the fields carry no energy density. In the general 
case, the decay rate depends on the equation of state in 
the background. 

Different behaviors of the spatial Stiickelberg fluctua- 
tion are related to the background temporal Stiickelberg 
field. This difference is the source of apparently conflict- 
ing claims in the literature regarding Stiickelberg per- 
turbation dynamics. The temporal Stiickelberg field in 
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the background is itself the choice of the time coordi- 
nate in which the massive gravity theory appears locally 
as the linearized Fierz-Pauli theory, i.e. it is a choice of 
unitary gauge for the covariant theory. For unitary time 
coinciding locally with isotropic time, the single propa- 
gating degree of freedom found in the decoupling limit 
[g, [l7| appears in the exact theory as well. For unitary 
time that scales with the spatial scale factor, there are 
no Stiickclbcrg dynamics, consistent with the open uni- 
verse solutions [201. Furthermore these solutions remain 
static not by possessing a vanishing stress tensor pertur- 
bation but rather by a potentially unstable balance of 
anisotropic stresses. 

The sign of the energy density carried by the 
Stiickelberg field ST, as well as the sign of the prefac- 
tor of the term in the quadratic action that generates 
the ST dynamics, are both determined by the quantity 
P[ , a constant determined by the parameters of the the- 
ory («3 and 0:4, see Eq. I30|) . P[ = —4 when = ay = 0, 
but can be either positive or negative for general 013 and 
ay. This is consistent with decoupling limit findings that 
perturbations around self-accelerating solutions generi- 
cally are ghost-like for 03 = 0:4 = 0, but can be made 
healthy in more general circumstances [(J [13] • P[ can be 
also be set to zero by a special choice. These models have 
no quadratic Stiickelberg contributions to the action and 
no linearized stress energy perturbations, consistent with 
the findings of Ref. [l,[lo[, who also made this parameter 
choice. 

In the exact theory, the dynamical system represented 
by the Stiickelberg fields exhibits some peculiar proper- 
ties that are obscured in decoupling limit analyses. The 
two Stiickelberg field perturbations require two sets of 
initial data - 6f(r,0) and 5T(r, 0) - and each obey cou- 
pled first order differential equations of motion. This 
is the amount of initial data and dynamics that we 
would expect for a single propagating degree of freedom. 
This situation might have been anticipated at the outset 
from a counting argument. Starting with 4 Stiickelberg 
scalar degrees of freedom, we remove two by restriction 



to spherical symmetry, and a third (the Boulware-Dcser 
ghost) is removed by construction in this theory. 

However, these two "half degrees of freedom are un- 
usual in their interrelation. As indicated above, the ST 
field has a first order equation of motion that is indepen- 
dent of the Sf field. The Sf dynamics do depend on the 
ST field, but they carry no stress energy at this order in 
perturbation theory. In addition to this, we find that the 
time derivatives of ST do not actually contribute energy 
or momentum density either, so that there is no obvi- 
ous classical instability associated with these fields even 
when they are apparently ghost-like. 

Finally, we uncovered a potential physical problem 
with certain background configurations for the temporal 
Stiickelberg field. The structure of the theory itself does 
not prevent the determinant of the fiducial metric from 
vanishing when solving the temporal or spatial boundary 
value problem. Some solutions for / pass through a zero 
determinant at finite radius (e.g. or time from well- 
posed initial values. At this point, the mapping between 
unitary and isotropic coordinates becomes singular. 

Since such a choice may have physical pathologies, we 
also constructed a new solution for / that passes from 
radiation domination through matter domination to self- 
acceleration without exhibiting this potential problem or 
resorting to static fields. We leave the larger question 
of the theoretical implications for the existence of poten- 
tially pathological solutions to a future work. 
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Appendix A: Second Order Action 

The quadratic Lagrangian for the perturbations (<5r, Sf, 8a, 8b) can be separated out as 

dsS'Mg + ^ST^g + 8\ T U g , 



8 z U a 



(Al) 



gU/ (Asm. 6) = a br ■ U/4 carries the relevant pieces of the metric determinant, S refers to pieces 



where U g 

involving Stiickelberg fluctuations and T those involving metric fluctuations. The Stuckclbcrg-Stuckclbcrg terms are 
1 



SssMg 



ST 2 ab [Pg + P['VX + P^'W ) + a 2 brPl,ST [5T'd g W + STd g W 



-a 2 brP[ST 5T'd g , + 8Td g + 8fd f + Sfd r VX + ^ar(2P 2 - x P! 2 )5f(STd r - 8T'd t ){ar) 



= P> 



-Sf8T eom + ydT 2 + a 2 brST(8T'd g , + STd g )VX 



(A2) 
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where we have used the background solution g = xgar, the polynomial recursion Eq. (|12[) . and the relationship between 
W and s/X on the background from Eq. (p~6|) . Note that equality here means equality up to total derivative terms. 
Here 

5T eom = d t {a 2 rb8TdjVx] + d r [a 2 rbSTd r Vx} - fmr[(ar)'ST - arST']. (A3) 
and its correspondence to Eq. (|20|) can be established by noting 



bVX 



f 
b 



df.Vx- 



i 



aVX 



r 



(A4) 



Next, the mixing between the Stiickelberg fields and the metric fluctuations is given by 



5 2 ST U g 



abr 



2P^ + P[{2 + ad a + ax d g >)VX + P^l + ad g >)W 



SaST 



+a 2 br 2 xo 



(P' 2 dgW + P[dgVX)Sa + {Pid 3 >W + P[d g ,VX)5a' 



-2a 2 br 2 P 2 Sa(STd g + 5T'd g >)W + a 2 r P^ + P[{1 + bd b )VX 



ST 
6b6T 



a 2 brP{{ {l + ad a )VX-x 



Sa 
a 



(A5) 



(Sa r)'d g , + Sa rd g (x VX -W)+ (1 + bd b )VX - x 



Both the 55 and ST pieces are directly proportional to P[. Note that 

f 2 - g 2 + b 2 W 



(1 + ad a )VX 



b 2 Vx 



(1 + bd b )VX 



- f 2 + a 2 W 
a 2 VX 



(A6) 



and these factors appear also in the perturbed energy and radial pressure relations, Eqs. (|21[) and (|23[) . since functional 
variation of the action with the metric must return these components. Note that there is no mixing between the 8f and 
metric fluctuations, as anticipated from our direct calculation of the stress energy tensor. On the other hand, variation 
of the action with respect to i5r yields an equation of motion for Sf that depends on all of the other perturbations. 
Finally, the metric-metric part is 

5 2 TT U g =3ar 2 P (bSa 2 + aSaSb) . (A7) 
As expected, it represents the perturbation to y/—g multiplying the effective background cosmological constant. 
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